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density of states of disordered semiconductor nanowires driven into a putative topological p-wave 
superconducting phase by spin-orbit coupling, Zeeman spin splitting and s-wave superconducting 
proximity effect induced by a nearby superconductor. Comparing with the corresponding theoretical 
self-consistent Born approximation (SCBA) results treating disorder effects, we comment on the 
topological phase diagram of the system in the presence of increasing disorder. Although disorder 
strongly suppresses the zero-bias peak (ZBP) associated with the Majorana zero mode, we find some 
clear remnant of a ZBP even when the topological gap has essentially vanished in the SCBA theory 
because of disorder. We explicitly compare effects of disorder on the numerical density of states in 
the topological and trivial phases. 
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I. INTRODUCTION 

The theoretical predictiorP^ that the combination of 
spin-orbit coupling, Zeeman spin splitting, and ordinary 
s-wave superconductivity could lead to an effective topo- 
logical superconducting phase under appropriate (and ex- 
perimentally achievable) conditions has led to an explo- 
sion of theoretical and experimental activities^ in semi- 
conductor nanowires (InSb or InAs) in proximity to a 
superconductor (NbTi or Al) in the presence of an ex- 
ternal magnetic field. The experimental finding^ of a 
ZBP, in precise agreement with the theoretical predic- 
tions in the differential tunneling conductance of an fnSb 
nanowire (in contact with a NbTiN superconducting sub- 
strate) at a finite external magnetic field (B ~ 0.1 — 1 T), 
followed by independent corroborative observatioiPHD] f 
such ZBP both in InSb and InAs nanowires in contact 
with superconducting Nb and Al by several groups, has 
created excitement in the condensed matter physics com- 
munity as well as the broader scientific community as per- 
haps the first direct evidence supporting the existence of 
the exotic, the elusive, and the emergent unparied Ma- 
jorana bound state in solids. Such excitement has in- 
variably been followed by a wave of skepticism as one 
would expect in a healthy and active scientific discipline 
with questions ranging all the way from whether such 
ZBP could arise from other (i.e. non-Majorana) origin 
to whether all aspects of the observed experimental phe- 
nomenology are consistent with the putative theoretical 
predictions on the topological superconductivity under- 
lying the existence of the Majorana mode. 

One particular issue, which is also the subject of 
the current work, attracting a great deal of theoret- 
ical attentiorP^HUl is the role of disorder in the Ma- 
jorana physics of superconductor-semiconductor hybrid 
structures. Disorder plays a key role in the Majo- 
rana physics because the underlying topological super- 



conducting phase hosting the Majorana mode (at de- 
fect sites) is essentially an effective spinless p-wave 
superconductor^^, which, unlike its s-wave counterpart, 
is not immune to non-magnetic elastic disorder (i.e. 
spin-independent momentum scattering) as was already 
known ten years ago 2 ^. Thus, even the simplest kind 
of disorder, namely zero-range random non-magnetic 
point elastic scatterers in the wire, could strongly af- 
fect the topological superconductivity in contrast to ordi- 
nary s-wave superconductivity which is immune to non- 
magnetic disorder, and the associated Majorana bound 
states by suppressing the (topological) superconducting 
gap^ and/or creating Andreev bound states in the super- 
conducting gap near zero-energy-H complicating the ob- 
servation and the interpretation of the ZBP. We do, how- 
ever, mention that elastic disorder or momentum scat- 
tering in the superconductor itself, no matter how strong 
(as long as it does not destroy the superconductor), does 
not affect the topological superconducting phase in the 
semicondutorJ^llln addition, it has recently been empha- 
sized that elastic disorder by itself could create a zero- 
bias peak (essentially, an anti-localization peak associ- 
ated with the disorder-induced quantum interference) in 
the non-topological phase in the presence of spin-orbit 
coupling and Zeeman splitting. Since the precise topolog- 
ical quantum critical point (as a function of the applied 
magnetic field) separating the topological and the triv- 
ial superconducting phase is in general not known in the 
experiments,^^ one cannot be absolutely sure that the 
observed ZBP is indeed a Majorana bound state (MBS) 
signature in the topological phase and not a trivial anti- 
localization peak in the non-topological phase. 

In the current work we consider disorder effects on 
MBS physics by directly calculating the density of states 
(DOS) of finite disordered nanowires in the presence 
of proximity-induced superconductivity taking into ac- 
count spin-orbit coupling and spin splitting arising re- 
spectively from Rashba and Zeeman effects in the wire. 
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We use random uncorrelated point scatterers with 5- 
function potential in the wire to represent the elastic dis- 
order. The theory follows the standard prescription^ of 
an exact diagonalization of the BdG equations in a mini- 
mal tight-binding model including superconducting pair- 
ing, Rashba spin-orbit coupling, and Zccman splitting in 
the Hamiltonian. The diagonalization of the discretized 
tight-binding Hamiltonian leads to the exact eigenstates 
of the system, which then immediately give the density of 
states (theoretical details are available in Ref. [26] and are 
not repeated here.) Comparison with theory in the pres- 
ence of disorder necccssitatcs the ensemble averaging over 
many different impurity configurations since each disor- 
der configuration produces its own unique result with 
random impurity-induced delta-function peaks in the su- 
perconducting gap. 



II. DOS IN THE DISORDERED TOPOLOGICAL 
PHASE 

In this section, we consider the effect of disorder on 
the SM/SC system starting from the topological phase 
in the clean limit. In Fig. 1 (with 8 panels, each repre- 
senting a different strength of disorder keeping all other 
parameters fixed), we show our numerical DOS as a func- 
tion of energy (E) for different relevant parameter sets in 
the topological phase of the system. We show both the 
ensemble averaged DOS using many-impurity configura- 
tions(but using the same disorder strength, i.e. the same 
impurity density and potential strength, changing only 
the random localitions of the impurities along the wire) 
and the typical DOS for a single impurity configuration 
in each case. Each panel corresponds to a specific dis- 
order strength (i.e. a fixed impurity density) and shows 
results for three different lengths of the nanowire. 

Since our calculations follow either Ref. 26 for the ex- 
act numerical treatment of Ref. |T3] for the SCBA theory, 
we refer the reader to those references for the technical 
details, which are actually pretty standardP^HUl 

The superconductor (SC)-semiconductor (SM) hybrid 
nanowire structure is characterized by a large number 
of independent parameters, both for the actual experi- 
mental laboratory systems and for our minimal theoret- 
ical model. The minimal set of parameters necessary to 
describe the system are the proximity-induced SC gap 
(A ) in the SM, the parent SC pairing potential (A s ), 
the Rashba spin-orbit coupling («_r), the spin splitting 
(Vz), the chemical potential (/z), the SM effective mass 
(to) which defines the SM tight-binding hopping param- 
eter, the nanowire length (L), and disorder (which we 
take to be the uncorrelated random white noise poten- 
tial associated with randomly located ^-function in real 
space spin- independent scattering centers). There are 
a few additional physically important parameters which 
are, however, not independent parameters of the the- 
ory: the coherence length in the nanowire (£), the SM 
mean-free path (A) due to disorder, and the number of 
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FIG. 1. Disorder averaged density of states in the topologi- 
cal phase for the semiconductor nanowire in a magnetic field 
with Zeeman splitting Vz — 5K, proximity- induced pairing 
potential of amplitude Ao = 3 K, Rashba spin-orbit cou- 
pling strength an = 0.3 eV — A. For this choice of param- 
eters the clean quasiparticle-gap A = 1.3 K and coherence 
length £ = 0.3 [im. The difference panels (a-h) correspond 
to different disorder strengths characterized by E s — h/2r = 
f 0, 50, fOO, 200, 460, 1300, 3300, 7300 mK. The corresponding 
mean-free paths are A are in the panels and the self-consistent 
Born gaps E SCB = 1.2, 1.1, 0.9, 0.8, 0.13, 0, 0, OK. 



occupied subbands (i.e. transverse quantized levels) in 
the nanowire which we take to be one throughout as- 
suming the system to be in the one-dimensional limit. 
(The single sub-band approximation, made entirely for 
the convenience of keeping the number of parameters in 
the model to be tractable, is a nonessential approxima- 
tion, and our qualitative results should be completely 
independent of this approximation.) In addition, there 
is an independent parameter defining the hopping am- 
plitude across the SC/SM interface which controls the 
proximity-induced SC pairing gap Ao in the SM in terms 
of the parent gap A s in the SC. Finally, the proximity 
gap (A) in the SM in the presence of spin-orbit cou- 
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pling and finite Zeeman splitting is reduced from Ao in a 
known manner. All the details for modeling the disorder 
are given in Ref. H3] where SCBA was used in contrast to 
our exact numerical diagonalization in the current work 
following Ref. One specific goal of our current work 
is to compare the analytic and simple SCBA theory^ 
with the exact tight-binding numerical analysis to test 
the limits of validity and the applicability of the SCBA 
theory which, being analytic, can be used rather easily. 

We choose parameters approximately consistent with 
the InSb/Nb systems studied experimentally irPl These 
are : Ao = 3K, an = 0.3eV — A (corresponding to an 
effective spin-orbit coupling strength m*a 2 R = 2.5K). 
Since out interest is in disorder effects, we focus on a 
range of magnetic fields with a spin splitting Vz ~ 5K 
(we vary it in a few cases only to change the proximity to 
the topological quantum phase transition point, separat- 
ing the topological and the trivial SC phase). Given that 
the condition for the topological SC phase to be realized 
in the SM is given bj3 Vz > A 2 + where A is the 
actual induced gap in the SM, the system should be in 
the MBS carrying topological phase for A < 5K (since 
fi = 0). Given that A < Ao = 3K, with the reduction of 
A below Ao arising from the existence of Vz ^ 0, our sys- 
tem is deep in the topological phase for the results shown 
in Fig.l since V z (= 5K) > A(= 1.3K) and n = 0. 

Each panel in Fig. 1 corresponds to a different disorder 
strength in the system, characterized by the correspond- 
ing level broadening E s — h/2r (where r is the scattering 
time — t — oo in the absence of disorder) or equivalently 
the mean free path A = vpr (where vf is the fermi ve- 
locity), both calculated in the SCBA according to Ref. 
13 for the given disorder in the wire. In each panel (and 
for each disorder) we show our DOS numerical results 
for two distinct wire lengths L = 1.5 fim and L — 3 \xm. 
In each case, we show both the ensemble averaged DOS 
results using an averaging over many (> 1000) random 
impurity configurations (keeping E s , A etc fixed) and the 
result for a typical single impurity configuration (the dis- 
tinct crosses or dots denoting delta functions for the DOS 
at the value of energy ). We emphasize that for an in- 
finitely long wire (L ^> £ s» 0.5 \im for our case) the DOS 
in the absence of disorder will vanish throughout the gap 
(±1.3 K in our case) with a (5-function peak at E = 
associated with the MBS at the wire edges. 

To characterize the disorder strength for the results in 
Fig. 1 (with panels (a) to (h) with increasing disorder 
keeping all other parameters fixed ), we use SCBA for 
this problem which was developed by us in Ref. 1131 The 
SCBA theory provides us with the SC gap in the topo- 
logical phase for a given disorder strength, allowing us 
to compare our direct (and exact) numerical calculation 
in the presence of disorder with the SCBA theory. We 
show the calculated SCBA gap in each case in the figure 
captions for the sake of direct comparison with the exact 
results in the figures. (The SCBA theory is obviously 
an ensemble averaged theory for the infinite system and 
does not depend on L.) It is clear from the result of Fig. 



1 that the analytical SCBA theory of Ref. H3] are in 
excellent qualitative agreement with the exact ensemble- 
averaged numerical results for the DOS even for E s s» A, 
where the topological gap essentially vanishes (panel (e) 
in Fig. 1) both according to the SCBA (i.e. E$cb = 0) 
and in our numerical results. While the disorder aver- 
aged DOS calculated by exact diagonalization does not 
stricly vanish inside the gap, the gap calculated with the 
SCBA can be identified with the peaks in the DOS at the 
edges of the gap. The closing of the gap within the SCBA 
coincides with the disappearance of the dips in the DOS 
around zero-energy. The DOS peak at E = associated 
with the MBS is continuously suppressed with increasing 
disorder, but quite amazingly there is a discernible DOS 
peak at E = Q even for E s (= 3.3K) > A(= 1.3K) where 
Escb = 0, and at best the topological superconductivity 
is gapless. 

The remarkable result, which is quite apparent in our 
Figs. 1 (e)-(g), is that the MBS peak of the DOS at 
E = is actually very robust to disorder and survives 
disorder strength substantially larger than that (typi- 
cally E s ~ A) destroying the induced superconducting 
gap A. Thus, in Figs. l(e)-(g), although the SCBA 
theory and our exact numerical results both show the 
system to be gapless with E s > A, the DOS peak at 
E = associated with the MBS persists until E s 3> A 
as in Fig. 1(h) where E s = 7X(> A = 1.3 K). It is 
not only that the MBS feature in the ensemble averaged 
DOS survives up to very strong disorder (e.g. the mean 
free path A = 0.5 /im, 0.2 fim and 0.06 \im respectively in 
Figs. l(e)-(g) which are smaller than the wire lengths 
of L = 1.5 /xm and 3/im used in our numerical work), 
the typical DOS for single random impurity configuration 
also shows peaks at E = as can clearly be seen in Figs. 
l(e)-(g) [ and as well as in Figs. l(a)-(d)], but not in Fig. 
1(h) where the very large disorder strength (A = 0.02 (j,m) 
suppresses both the ensemble averaged MBS peak as well 
as the single configuration peak at E = 0. The survival 
of the zero-energy DOS peak well above the point where 
the SC gap is completely suppressed by disorder is an im- 
portant new result of our exact numerical work directly 
establishing the possible theoretical existence of a gapless 
topological SC phase. 



A. Griffiths effects 

The origin of the E = peak in the DOS in the 
strongly disordered case can be related to the Grif- 
fiths effect previously considered for the spinless p-wave 
superconductor^. The Griffiths effect in the semiconduc- 
tor nanowire at finite Vz arises from the disorder-induced 
variation of the chemical potential, which can lead to a 
transition from a topological phase to a non-topological 
phase and vice- versa. The variation of the effective chem- 
ical potential can lead to domain-walls between topologi- 
cal and non-topological regions each of which would sup- 
port a local zero-energy MBS. Since each region is of 
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finite extent, the MBSs are a finite distance apart and 
split into conventional states with a non-zero energy. In 
fact, by carefully considering the distribution of the dis- 
tances between the MBSs, it has been showrP^l that this 
splitting typically leads to a singular peak in the DOS at 
E = 0, which diverges as E — > as a power-law in E. 
The topological phase is of course characterized by true 
zero-energy edge states, whose energy is exponentially 
small in the length of the system L. Thus, in the long 
wire limit, the topological phase and the non-topological 
phase both have power-law divergent DOS due to the 
Griffiths effects, but the topological phase has a pair of 
zero-modes exponentially close to zero energjG^. This 
distinction between exponential versus power-law in the 
length seems to appear in the DOS plotted in Fig 1. At 
small disorder, there is a sharp peak which is very weakly 
split. At larger disorder (panels (e,f,g,h)) , which is ap- 
proximately when the gap vanishes i.e. E$cba 0, the 
E ~ peak in the DOS is broadened. 

For finite length systems, the sharp transition between 
the topological and non-topological phase at finite disor- 
der becomes a crossover. While the topological invariant 
can be co mputed e ven for a disordered, but strictly infi- 
nite system ^ 23 ! 24 !, we have not done so in the present 
work because we restrict ourselves to systems of wire 
lengths comparable to the experimental systems. For 
such finite wires the distinction between the topological 
and non-topological phase is not sharp. In fact, the zero- 
modes arising from the Griffiths effect are themselves 
MBSs. The MBSs characterizing the topological phase 
are special only in the sense that they occur near the 
ends of a finite system and are therefore separated from 
other low-energy MBSs by a distance of the order of the 
length of the wire. This can also occur from the Grif- 
fiths effect in some realizations of disorder, even in the 
non-topological parameter regime. Therefore, the Grif- 
fiths singularities seen in the disordered Fig. 1 indicate 
the presence of several low-lying MBSs in the spectrum 
of the finite wire. For very large disorder the system be- 
comes completely non-topological (and gapless) as in Fig. 
1(h). 

There are three important messages following from our 
DOS results for the topological phase in the presence of 
disorder: (1) the zero-energy DOS peak associated with 
the MBS is strongly suppressed by disorder; (2) SCBA 
is an excellent quantitative approximation for calculat- 
ing the SC gap in the topological phase including effects 
of disorder; and (3) most importantly, the zero-enenrgy 
MBS peak is very robust against disorder and survives 
well after the SC gap in the topological phase has been 
suppressed to zero, and disappears only when the mean 
free path A < £/6 or E s > 6 A (for the systems are stud- 
ied) whereas the topological SC gap vanishes for A < £. 

Given the quantitative validity of the SCBA, we can 
calculate the phase diagram of the SC/SM structure in 
the presence of disorder by using the analytical SCBA 
theory.^ We show our SCBA-calculated quantum phase 
diagram of the system in Fig. 2. In Fig. 2, we fix all pa- 
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FIG. 2. Quasiparticle gap A (in the color bar) versus Zeeman 
potential Vz and scattering rate E a in the topological phase. 
The black region represents the gapless phase and is therefore 
not topologically robust. 



rameters except for disorder (E s ) and spin-splitting (Vz) 
with the color representing the calculated SC gap E$cba 
in the presence of disorder (by definition Escba = A 
for E s = 0). The black region represents gapless super- 
conducivity. We note, based on our numerical results 
of Fig. 1, that a large part of the black region in the 
phase diagram allows for well-defined zero-energy DOS 
MBS peaks although the system is essentially a gapless 
topological superconductor in this regime. Much of this 
SCBA gapless regime is dominated by Griffiths physics 
except for very large disorder when the system eventu- 
ally becomes nontopological (i.e. even the finite topo- 
logical segments disappear). Much of this SCBA gap- 
less regime is dominated by Griffiths physics except for 
very large disorder when the system eventually becomes 
non-topological. Except for our current results showing 
the well-defined robust persistence of the MBS peak even 
in the gapless topological regime (the black region in 
Fig. 2), one would have concluded that SCBA predicts a 
rather gloomy picture for the existence of the Majorana 
mode in the disordered SC/SM hybrid structure since, 
without our current exact results, the conclusion would 
have been that there cannot be any MBS in the black 
region of the phase diagram. Of course the quantitative 
details of the SCBA phase diagram depend on the SO- 
coupling strength and the topological phase with a large 
SC gap is easily achieved by large (small) SO coupling 
(disorder). 



III. DOS IN THE DISORDERED 
NON-TOPOLOGICAL PHASE 

Next we comment on the effect of dis order in t he non- 
topological SC regime, i.e., for Vz < \J & 2 + f° r the 
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FIG. 3. Disorder averaged density of states in the non- 
topological phase for the semiconductor nanowire in a mag- 
netic field with Zeeman splitting Vz — 0, 1 K (panels (a, b) 
respectively) at n — for different disorder strengths charac- 
terized by scattering rates E s = h/r. 



FIG. 4. Disorder averaged density of states in the non- 
topological phase for the semiconductor nanowire in a mag- 
netic field with Zeeman splitting Vz = 4, 5K (panels (a, 6) 
respectively) at /i = 5 K for different disorder strengths char- 
acterized by scattering rates E s — h/r. 



SC/SM hybrid structure. In Fig. 3(a) and (b), we show 
our calculated numerical DOS in the non-topological SC 
phase (for fj, = 0) for Vz = (3a) and IK (3b) in the 
presence of disorder. All parameters other than Vz are 
exactly the same as in Fig. 1. We see that, at least 
for /i = 0, the behavior of the DOS is similar to a clas- 
sic s-wave SC (even for Vz = 1 K which simply reduces 
the gap from Ao = 3 K to 2K) with essentially no dis- 
cernible effect on the DOS. Even for disorder as large as 
E s > 7 K > 2Ao, we do not see any structure developing 
in the SC DOS gap which remains completely flat. For 
Vz = (Fig. (3a)) this is a direct consequence of Ander- 
son's theorem where the robustness of the gap arises from 
time-reversal symmetry. We see in Fig. (3b) that this be- 
havior persists for small Zeeman fields as long as Vz A. 
This behavior can be expected based on previous studies 
on the bound states of single impurities in spin-orbit cou- 
pled nanowires in proximity to superconductors^. There 



it was found that the low-energy sub-gap states appear 
for short-ranged impurities only in the topological phase. 
However, this conclusion might not apply to longer-range 
disorder because in principle, at any non-zero value of 
Zeeman potential puddles can lead to the formation of a 
pair of low-energy MBSs. In our numerical work here we 
restrict ourselves only to zero-range white noise disorder 
in the semiconductor. 

The situation, however, changes qualitatively when we 
consider the limit Vz ^ A. We show these results in 
Fig. 4 where V z — 1 K, fi = 5 K and V z = /i = 5 K 
are shown for several values of the disorder parameter 
E s (= H/2t). All other parameters are the same as in Fig. 
1 and 2. We note that the situation in Fig. 4 describes 
the non-topological phase since Vz = A < \/^ 2 + A 2 . 
In both panels, large disorder has a strong effect on the 
SC DOS shrinking the SC gap considerably. However, 
in Fig. 4(b), where Vz > A, eventually for E s = 0.2 K 



6 




FIG. 5. Quasiparticle gap (A) versus Zeeman potential Vz 
for different values (fi — 0, 5K) for the chemical potential [i. 
The vanishing of the quasiparticle gap marks the topological 
quantum phase transition from the non-topological phase at 
small Vz to the topological phase at large Vz- 



and 1.3 K (and V z = A* = 5K), the DOS seems to "flip" 
and the dip at E = becomes a peak at E = with 
a concomitant vanishing of any SC gap feature in the 
data. For even larger E s (E s = 7 K in Fig. 4) the peak 
feature in the DOS is suppressed, but there is a clear 
E = DOS peak for intermediate (but still large with 
E s > A) disorder (E a = 0.2 K and 1.3 K in Fig. 4) 
where the SC gap has disappeared, but a peak has de- 
veloped in the DOS at zero-energy. We believe that the 
zero-energy DOS peak in Fig. 4(b) for E s = 0.2 K and 
1.3 if has the same origin as the physics recently dis- 
cussed in several publications^^. The hallmark of this 
"trivial" DOS peak (arising from the competition among 
spin splitting, spin-orbit coupling, and superconductiv- 
ity) are that (1) it arises only for large Zeeman splitting 
in the non-topological phase; (2) it occurs only after the 
SC quasiparticle gap has been completely suppressed by 
disorder; (3) it exists only in an intermediate disorder 
range where the superconducting quasiparticle gap has 
been suppressed, vanishing for larger disorder and the 
peak becoming a dip (i.e. the SC gap) for smaller disor- 
der. 

In order to specify where precisely in the phase dia- 
gram we are obtaining our exact numerical results in the 
presence of disorder we finally show in Fig. 5 our cal- 
culated SC gap as a function of Vz (for zero disorder) 
for the parameters chosen in our calculations. For fi = 
(i.e. Figs. 1-3), the topological quantum phase transition 
(TQPT) happens at Vz = A (= 3K in our choice), and 
our Fig. 1 belongs to the topological phase (Vz = 5K, 
to the right of the TQPT) whereas our Fig. 2 belongs 
to the non-topological (Vz = 0,1 K) to the left of the 
TQPT. Results for Fig. 4 are obtained for Vz = A* (and 
Vz < £t) which never can manifest a TQPT since the 



Vz = \j [i 2 + A 2 condition cannot be satisfied except for 
A = (i.e. Vz —> oo)- but, even the nontopological SC is 
strongly affected by disorder (for Vz 3> A) here since the 
combination of spin-orbit coupling and Zeeman splitting 
makes the Anderson theorem moot. Thus the DOS peak 
in Fig. 4, associated with anti-localizationf^HUl arises 
purely in the trivial phase with a completely suppressed 
quasiparticle gap which might enable its experimental 
distinction from the MBS peak in the topological SC 
phase. 



IV. DISCUSSION 

Finally we mention that the spin-orbit coupled semi- 
conductor system thus has two distinct topological quan- 
tum phase transitions in the presence of the supercon- 
ducting proximity effect and Zeeman splitting. The first 
one is driven by the Zeeman field as originally predicted 
by Sau et aPwith the TQPT defined by V z = \J + A 2 
assuming a low disorder situation. The second one is 
driven by increasing disorder (E s ) in the finite Zee- 
man splitting situation (i.e. V z > \/a* 2 + A 2 ) where 
the topological SC phase is destroyed by disorder (for 
E s > A), leading to a gapless non-topological phase dom- 
inated by the Griffiths physics as originally envisioned by 
Motrunich et aP^l Since the effective SC gap for Vz ^ 
depends on Vz, and in particular, A cx V^ 1 for V z A 
[see Refs. l5l and [TBI for details], we expect that there are 
two distinct magnetic field driven TQPTs in the semicon- 
ductor nanowires - - the first one is the TQPT predicted 
in Rcf. [2HS1 taking the system from the trivial s-wave SC 
(induced by proximity effect) to the (effective p-wave) 
topological Majorana carrying p— wave SC phase at low 
disorder (i.e. A 3> E s ), and then the second one at much 
higher Zeeman splitting (so that A < E s ) where disorder 
drives the system from a gapless topological SC to a non- 
topological Griffiths phase. It is unlikely that this second 
(purely disorder driven) TQPT would be experimentally 
accessibly since the gapless nature of the SC phase (i.e. 
black region in Fig. 2 or the situation corresponding 
to Fig. 1 (e)-(h)) would make the finite temperature of 
the experimental system behave like a very high tem- 
perature (i.e. T> A), making any experimental study 
of this disorder-driven TQPT difficult, if not impossible. 
Thus, the effective gapless nature of the system in Fig. 
l(e)-(g) would make it very unlikely that the DOS peak 
(which is quite obvious in our theoretical results) could 
be studied experimentally. The best hope for the direct 
experimental study of the MBS physics is therefore to 
have a large SC gap (as in Fig. l(a)-(d) or in the non- 
black region of Fig. 2) in the topological phase which 
necessitates having low effective disorder (E s <C A), a 
condition guaranteed by having very clean semiconduc- 
tor wires and or/very strong SO coupling in the material. 
We add here that our theoretical results presented in this 
paper cannot be compared (or connected) with the ex- 
periments in any way since we only calculate bulk DOS 
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of the system, which cannot be directly probed experi- 
mentally. 

V. CONCLUSION 

In summary, we have studied the disorder averaged 
DOS of a disordered spin-orbit coupled nanowire in prox- 
imity to a superconductor in both the topological and 
non-topological parameter regime. The features in the 
DOS associated with the superconducting gap in the 
topological phase appear to be in good quantitative 
agreement with the SCBA from previous work^l. Consis- 
tent with previous results^, we find that the dips in the 
DOS associated with the quasiparticle gap disappear for 



relatively modest (E s > A) amounts of disorder. How- 
ever, a zero-energy peak associated with MBSs generated 
by the Griffiths effect survives to much higher levels of 
disorder. The DOS peak associated with different levels 
of disorder arising from the Griffiths effect starting from 
the topological phase and the antilocalization peak start- 
ing from the non-topological phase appear to be different 
enough that one might distinguish them qualitatively. Of 
course, at this point one does not expect a sharp distinc- 
tion between the peak arising from the Griffith's effect 
and the antilocalization effect because they are zero en- 
ergy peaks in the DOS in the non-topological phase in 
the same symmetry class i.e. class D. 

This work is supported by Microsoft Q, JQI-NSF-PFC, 
and the Harvard Quantum Optics Center. 
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